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Abstract. - The possible "phase diagrams" for shear-induced phase transitions between two 
phases are collected. We consider shear-thickening and shear-thinning fluids, under conditions of 
both common strain rate and common stress in the two phases, and present the four fundamental 
shear stress vs. strain-rate curves and discuss their concentration dependence. We outline how 
to construct more complicated phase diagrams, discuss in which class various experimental 
systems fall, and sketch how to reconstruct the phase diagrams from rheological measurements. 



Introduction. - There is growing interest in the flow behavior of complex fluids, including 
worm-like micellar surfactant solutions jl|-|8| , lamellar and "onion" surfactant systems |)] |ll| , 
and liquid crystals [[l2| |l3| . In these and other systems the microstructure is altered by 
flow such that a bulk transition, reminiscent of an equilibrium phase transition, can occur. 
Signatures of this transition are kinks, plateaus, or non-monotonic behavior in the measured 
"flow curve" , that is, the relationship between applied shear stress and mean strain rate; 
structural changes observable by X-ray ]|, [l2|], neutron or light scattering; explicit 

measurements of macroscopically discontinuous flow profiles ||; and visual confirmation of 
phase separation and coexistence |||7j 



Despite the surfeit of experiments, theories have been limited to a few systems (micelles 14 
[15| and liquid crystals (|l6|— |l9|] ) , and few have attempted to calculate an entire phase diagram 
for a complex fluids solution in flow p9| . In a remarkable paper, Schmitt et al. classified the 
instabilities that can occur in complex fluid solutions, and clarified the relationship between the 
nature of the instability (either in primarily the concentration, "spinodal" , or the strain rate, 
"mechanical" ) , the shape of the flow curves, and the orientation of the interface that initially 
develops during the instability |6|. However, as noted in the orientation of an interface 
in an initial instability may or may not be relevant to the orientation of the macroscopic 
phase separated state. In this work I outline the different macroscopic phase diagrams that 
can occur in complex fluid solutions in planar shear flow and describe how phase diagrams 

(*) p .d. olmstedOleeds .ac.uk 

Typeset using EURO-TfX 



2 



EUROPHYSICS LETTERS 



determine "flow curves" , the relation between applied shear stress and measured mean strain 
rate. No attempt is made to pose or solve any specific dynamical models (see ^0|), but 
rather to explore the consequences of possible phase diagrams and provide a phcnomcnological 
framework within which to understand the rapidly growing body of rheological experiments. 



Phenomenology of Phase Coexistence in Shear Flow. - Calculations of phase behavior 
requires first determining the relevant microstructure of the quiescent and flow-induced phases, 
and deriving equations of motion for the appropriate variables (momentum, mass concentration 
(j), and variables such as liquid crystalline or crystalline order, or aggregate shape and size 
distributions). These equations of motion determine the steady state flow curves corresponding 
to different microstructures. For example, fig. [l] shows a hypothetical shear-thinning fluid I 
which can be sheared into either a less viscous fluid Si, or a gel-like state S2- For this particular 
fluid the shear-induced phases are stabilized only in finite flow; systems near equilibrium 
phase transitions such as liquid crystals may have multiple stable flow curves in the limit of 
zero flow while other systems such as dilute wormlike micelles [Q and "onions" |)-1I 
apparently have flow branches which are only stabilized in finite flow. The task of theory, and 
indeed of experiments, is to understand how material on different flow curves may coexist, and 
what controls the stability of one phase with respect to another. 

At coexistence a fluid partitions into phases with, in principle, different concentrations. In 
the fluid of fig. [I] the I and Si states could coexist at the same shear stress (which would require 
the interface between phases to lie in the velocity- vorticity plane); as could the Si and S2 states; 
the I and Si phases could not coexist at a common shear stress. Phase coexistence among all 
three possible pairs of states is conceivable at a common strain rate (requiring an interface lying 
in the velocity- velocity-gradient plane). Phase diagrams may be constructed by determining 
the states on distinct stable flow curves for which the driving force for particle exchange 
vanishes, and for which a mechanically stable stationary interfacial solution between phases 
can be found. This procedure has been developed for dynamical models which incorporate both 
smooth |^,[l8|,[l9| and sharp |^0|,|l| interfaces. In the former case gradient terms are included 
in the equations of motion, while in the latter case a mechanical condition on interface stability 
is required. Complete phase diagrams have been calculated for a model system of rigid rod 




Fig. 1. - Steady state constitutive curves 17(7, (j>) for a hypothetical fluid with a shear-thinning qui- 
escent state I, a shear-induced phase Si, and a shear- induced "gel" S2 with an apparent yield stress. 
All flow branches have concentration dependences as shown, and are only stable at the points shown. 
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Fig. 2. - Fundamental rheological signatures for phase separation at common stress (top row, Al and 
A2) or common strain rate (bottom row, Bl and B2), for essentially shear thinning (middle, Al and 
Bl) or shear thickening (right, A2 and B2) transitions, a and 7 are the imposed (measured) shear 
stress and strain rate, respectively. The thick solid lines denote the two phase region if 8<f> — 0, and 
the broken lines denote the two-phase region in the event of a finite width biphasic region \8<f>\ > 0. 

suspensions, for both common stress and common strain rate geometries ||l9| . Unfortunately, 
we still lack methods for choosing from several candidate possibilities for phase coexistence. 

To begin, consider a complex fluid solution which possesses two distinct phases under shear 
flow: we denote these I and S (for the example fluid of fig. |l| these could be any of the I, Si 
and S2 states). We assume this fluid has steady-state macroscopic phase coexistence observed 
experimentally and hence has a phase diagram which may, in principle, be calculated from the 
relevant dynamical equations of motion. For simplicity we ignore interesting complications 
associated with secondary-flow or dynamical instabilities which can induce non-stationary 
oscillating or chaotic steady states; and ignore any effects of curved geometries. 

Given such a fluid, there are two possible geometries for phase coexistence: (A) common 
stress phase separation, for which phases coexist in the flow gradient direction, and (B) 
common strain rate phase separation, for which phases coexist in the vorticity direction. For 
each case, flow can induce a transition to cither a less viscous phase (Al, Bl) or a more 
viscous phase (A2, B2). These geometries and flow curves are shown in fig. |[ To understand 
the concentration dependence of these flow curves we must examine the entire phase diagrams: 
some possibilities are shown in fig. I Upon increasing the concentration the S phase may be 
stabilized at either weaker flows (Al, Bl, A2, etc.) or stronger flows (Al', Bl', A2', etc.). 

Common Stress Phase Separation. — We first consider coexistence at a common stress (Al, 
A2, Al', A2'), for which the mean strain rate 7 and mean concentration <p at stress a are 
partitioned according to 



7 = aji{^i) + (1 - a) 7s OAs), 
<(> = a (}>r + (1- a) <f>s, 



(1) 

(2) 
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with a the fraction of material in the I phase. Phase coexistence occurs in a region in the a—<f> 
plane, with horizontal tie lines connecting the coexisting points (0/, jj) and (0s, is)- This may 
also be represented in the 7— <f> plane. The lines {(f>i(a), 4>s(cr),'yi(^)i(a)),'ys(4's(o'))} bound the 
phase coexistence regions. The slope of tie lines in the 7— plane reflects the compositions of 
the two phases; vertical tie lines imply equal concentrations, while a non-infinite slope implies 
different concentrations. The flow curves (7(7,0) can be calculated from the phase diagrams 
and eqs. |l| and |^ by applying the lever rule. They are non-analytic at the stresses and strain 
rates which bound the biphasic region, and would be obtained in experiment at a given mean 
concentration 0. In the biphasic region the stress increases by crossing successive tie lines in 
the 7 — plane, and varies according to the tie line spacing and splay [pj| . From eqs. [l] and || 
one can calculate the shape of the stress "plateau" B: 



da 



a 

m 



Vs 



m(<7, 0) 



1 



a 



Ts'ns 7/?7/ 



(3) 



where m(a, 0) is the slope of the tie line, 77/. = da/dj is the local viscosity of the kth branch, 
and "/' k = d^k/dc/). Experiments on worm-like micelles have revealed that some systems can 
be "superstressed" to a metastable state, with the stress entering the biphasic region |^,^,^|. 




Fig. 3. - "Fundamental" phase diagrams and flow curves. These have not been calculated, but would 
be calculated from the homogeneous flow curves (e.g. fig. |l|) by the methods of fij| or @||. The 
phase boundaries 75, ji, as, 07 in the top two rows are coincident with the kinks in the flow curves. 
Top row: shear diagrams in the 7— <f> plane. Middle row: shear diagrams in the a— <f> plane. Thick lines 
in both planes correspond to the I phase boundary. Bottom row: flow curves (7(7) for representative 
mean concentrations (j>i and 02, with thick lines denoting the phase coexistence region; in each case 
0i has a narrower window of coexistence in the relevant control variable plane than does 02. Tie lines 
in the 7 — plane move to higher strain rates with increasing stress. Al and A2: phase separation 
at common stress with the S phase stabilized at a lower stress for increasing concentration. Al' and 
A2': phase separation at a common stress, with S phase stabilized at a higher stress for increasing 
concentration. Bl and B2 phase separation at a common strain rate. 
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For an Al (or AT) shear-thinning transition the strain rate increases as the stress is 
increased through the two phase region. The stress through the two-phase region increases 
less for mean concentrations (f> f° r which the biphasic window is narrow {e.g. fa) or the tie 
lines arc steep in the 7—0 plancl, and deviates more from constant for a wider biphasic region 
{e.g. (f> = fa) with more vertical tie lines. Hence, a signature of substantial concentration 
difference 5(f> = <f>j — (f>s between coexisting phases is a large increase in stress through the 
biphasic region. This is consistent with calculations on a model for rigid rods in shear flow [jl9| , 
and measurements on shear-thinning wormlike micelles . In dilute micelles 5(f> is typically 
negligible and the stress plateau is nearly flat; while 5(f> increases for concentrated solutions 
near an underlying isotropic-nematic transition and the stress plateau acquires a shape. 

For the AI fluid the S phase can have a larger (m(a, fa > 0) or smaller (m(<7, fa < 0) strain 
rate than the I phase. In the latter case the S phase has a larger effective viscosity, and for 
a small enough 5(f) the fluid crosses over to a characteristic shear-thickening fluid, A2. In this 
case the strain rate actually decreases as the stress is increased through the two phase region. 
That is, the system remains in the I phase below a critical strain rate at which a band of 
more viscous material S develops whose small strain rate reduces the measured strain rate. 
Hence the system enters the biphasic region in the 7 — cj> plane at the high strain rate, and 
traverses from top to bottom. Upon increasing the stress the system converts into the S phase 
with, in general, changing I and S concentrations, until the I phase disappears at the lower 
boundary of the two-phase envelope in the 7 — <f> plane. For higher stresses the system follows 
the constitutive branch for the S phase and traces a vertical path in the 7 — plane. Hence 
the flow curve (7(7) has a distinct S shape. As with flow Al, the shape of the flow curve in 
the biphasic region depends on the slope of the tie lines {i.e. 5fa in the 7 — </> plane: vertical 
tie lines 5(f) = imply a vertical jump in 17(7), while finite 5(f) and flatter tic lines imply a 
gentler slope for (7(7). Flow A2' has been seen in shear-thickening wormlike micelles 0, and 
calculated using a phenomenological model |po|] . Flow A2 (or A2') has been seen in surfactant 
onion crystals under shear, although the band geometry has not been verified fl(i|| . 

We have drawn (Al,A2,Al',A2') with finite biphasic regions at zero stress, corresponding 
to a perturbation of an existing phase transition. However, fluids such as the example 
fluid in fig. [l] would have phase coexistence only above a finite stress. The construction of 
these phase diagrams implies that tie lines move to higher strain rates with increasing stress 
("conventional"). Although intuitive, the reverse ("unconventional") does not violate any 
physical laws. For example, a version of A2' with unconventional positive-sloped tie lines, 
would imply a a—cf> biphasic region moving to smaller (f> with increasing a. However, a version 
of A2' with negative-sloped tie lines yields a nonsensical phase diagram. 

Common Strain Rate Phase Separation. - For coexistence at a common strain rate the 
shear stress is partitioned according to 



and the two main classes of flows are shown as Bl and B2. We do not show analogous diagrams 
Bl' and B2', for which the coexistence strain rate increases with increasing concentration. For 
flow B2 the stress increases through the two phase region and the transition is to either a 
thicker or (for large enough 5(f), as with flow Al) thinner phase. The shape of the stress 
through the two phase region is given, from eqs. |^ and |] 



a = a a 1 (fa) + (1 - a) a s ((f>s), 



(4) 



da 

(l A ; 



rji a + T)s (1 - a) - 771(7, fa 




(5) 



where m(j, fa is the slope of the tie line and a' k = d<Jk/d(f>. In the limit of no concentration 
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difference (5(f) = or 771(7, 4>) — <J ( A i) is vertical through the two phase region. Flow Bf is 
a shear thinning transition, and the stress decreases through the two-phase region. Banding 



geometries and flow curves similar to Bl were reported in an onion system 
may not have been steady state. B2 has been reported in an onion system 



although this 



Applications. - The shear diagrams presented do not exhaust the possibilities; these flow 
behaviors can be combined smoothly in many ways. We present two possibilities in fig. ^| 
Fluid Cf phase separates at a common stress, and 5<p widens with increasing concentration, 
crossing over from shear thickening (A2) to shear thinning (Al) as the concentration increases. 
Fluid C2 phase separates at a common strain rate, with a coexistence width S(f> that narrows 
with increasing </>, combining Bf and B2' behavior. 

Typical experiments extract a series of steady state flow curves for different concentrations 
{4>i}. Consider coexistence at common stress. The flow curves have kinks at the boundaries of 
the biphasic region, {(era, jn, 4>i)} and {(csi, Tsi, 4>i)}, which should be corroborated optically 
or otherwise. The strain rates and concentrations of the coexisting states may be determined 
as follows. One first fits curves <t/(<^), (ts(4>), 7/(0), 7s (</>) to the values extracted from the 
kinks. A horizontal tie line connecting (</>/, cf>s) at a given stress <r* may then be read off the 
biphasic boundaries, 01 (</>/*) = us(4>s*) = in the a— 4> plane; and the corresponding tie lines 
on the 7 — </> plane connects the points 7/ (</>/*) and js((f>S*)- In this way the complete shear 
diagrams and characteristics of coexisting states, may be constructed. As a check, the shape 
of 17(7) through the two-phase region may be computed by traversing the biphasic regions of 
the and a— <fi planes and using the lever rule to construct the mean strain rate via eq. ([!]). 
This analysis is necessary if the coexisting concentrations cannot be determined directly and, 
in the two phase region, the flow curve has a slope appreciably different from zero (common 
stress phase separation) or infinity (common strain rate phase separation). 

Returning to the hypothetical fluid of fig. |l|, we expect common stress phase separations 
I-Si and S1-S2 to be classes Al and A2 (or Al' and A2'), respectively, with the former 




Fig. 4. - More complex possibilities for phase separation at common stress (CI: top row) or common 
strain rate (C2: bottom row). <pi,cj>2, and 03 are representative mean concentrations. 
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disappearing at higher stresses; and common strain rate phase separations I-Sj., I-S2, and 
S1-S2 to be classes Bl, B2, and B2 (or Bl', B2', B2'), respectively. Since flow is necessary to 
stabilize the Si and S2 phases, the biphasic window in this system would appear at a finite 
stress or strain rate. Note that the accessibility and even stability of composite flow curves 
depends on the control variable of the rheometer. For example, composite flow curves with 
negative slope dcr/dj < (e.g. A2 and Bl) could be mechanically unstable ||, although ref. 
accessed such a curve under controlled stress conditions. Finally, kinetic possibilities such as 
metastability and hysteresis are sure to enrich the relatively simple scenarios proposed above. 

It is a pleasure to thank David Lu, Armand Ajdari, Ovidiu Radulescu, Jacqueline Goveas, 
and Tom McLeish for encouragement and enjoyable discussions. 
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